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Biphoton compression in standard optical fiber: exact numerical calculation 
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Generation of two-photon wavepackets, produced by spontaneous parametric down conversion in 
crystals with hnearly chirped quasi-phase matching grating, is analyzed. Although being spectrally 
broad, two-photon wavepackets produced this way are not Fourier transform limited. In the paper 
we discuss the temporal compression of the wavepackets, exploiting the insertion of a standard 
optical fiber in the path of one of the two photons. The effect is analyzed by means of full numerical 
calculation and the exact dispersion dependencies in both the crystal and the fiber are considered. 
The study opens the way to the practical realization of this idea. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

One of the problems in modern quantum optics is gen- 
eration of two-photon light with given spectral and tem- 
poral properties. In particular, of considerable interest is 
preparation of two-photon wavepackets with small corre- 
lation times. Studies in this direction are important, first 
of all, from the viewpoint of the fundamental question: 
how well can we localize a photon of a correlated pair 
if its twin is registered at a certain time? Furthermore, 
two-photon wavepackets with small correlation times are 
required for quantum technology applications involving 
interactions between matter and nonclassical light UJ 
[Bi], as well as for quantum metrology applications 
Several ideas have been suggested for generating two- 
photon states with broad spectra, all based on Spon- 
taneous Parametric Down Conversion (SPDC). Among 
them, one can mention prisms or diffraction gratings in- 
troducing a frequency chir p M , SPDC in aperiodically 
poled crystals (APPC) jsl-llQj or in crystals with tem- 
perature gradients However, a broad spectrum of 
two-photon light does not necessarily lead to small cor- 
relation times, although the inverse is true p!2l-[T^. Sim- 
ilarly, in classical optics, a broadband pulse is not always 
short in time, although the spectrum of a short pulse is 
always broad. Indeed, the spectrum broadening intro- 
duced in Refs. [3, |lO|, lU is inhomogeneous, in the sense 
that two-photon wavepackets generated in different parts 
of the crystal have different spectra. As a result, the 
two-photon spectral amplitude has a phase (a frequency 
chirp) depending nonlinearly on the frequency and mak- 
ing the two-photon wavepackets not Fourier transform- 
limited. Therefore, in these examples the biphotons are 
not short in time despite their broad frequency spec- 
trum. In Ref. [9] it was mentioned that time compression 
of such two-photon wavepackets requires compensation 
for their frequency chirp, although the way to eliminate 
the chirp was not specified. Recently we have demon- 
strated [15] that broadband biphoton wave packets pro- 
duced via SPDC in APPC can be compressed in time by 
using group-velocity dispersion (GVD) in standard opti- 



cal fiber. A Fourier limited compression was obtained in 
[m via an analytical calculation within the approxima- 
tion of large difference in the group velocities of the two 
SPDC photons and neglecting third-order GVD of the 
fiber. Furthermore, we have given an example of a com- 
plete numerical calculation demonstrating that higher or- 
ders do not change this result. 

With the purpose of paving the way for a realistic ex- 
perimental implementation of this idea and application 
to quantum technologies [16], here we present and discuss 
a systematic study of this phenomenon. Furthermore, we 
present some optimal situation where the full numerical 
calculation shows an effect that can be clearly observed 
with a realistic set-up. 



II. SPECTRAL PROPERTIES OF BIPHOTONS 

EMITTED IN A MEDIUM WITH 
SECOND-ORDER SUSCEPTIBILITY VARYING 
ALONG THE PUMP PROPAGATION 

Spontaneous Parametric Down Conversion occurs in 
non-linear dielectric crystals and consists of the spon- 
taneous decay of a photon (pump) of frequency ujp in 
two photons of lower frequencies uJs and cj^, historically 
named signal and idler, satisfying, in an ideal case, the 
constraints of energy and momentum conservation. In 
what follows we will consider type II collinear SPDC in- 
teraction, with the pump being a linearly (ordinarily) po- 
larized monochromatic plane wave propagating along the 
z direction and the signal and idler fields being extraordi- 
narily and ordinarily polarized, respectively. Considering 
signal and idler fields initially in the vacuum state, the 
output state is 

\^) = |0) + j du;sduj^F{ujs,uJi)alhl\0), (1) 

where, due to the weak interaction, only the first two 
terms of the perturbative expansion are considered pjf . 
In Eq. ([1]), aj^^ and are the photon creation operators 
for signal and idler photons, respectively. 
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The function F{(Js^^i) is the Two-Photon Spectral 
AmpHtude (TPSA) and accounts, in particular, for the 
non-linearity of the crystal x*^^^ and the wavevector mis- 
match kp — kg — ki. The second-order susceptibility is 
assumed to depend on the pump propagation direction z 
only; for cw pumping, the TPSA is 



(2) 

where L is the length of the crystal and the wavevectors 
kj = ^^^^^^^^ (j = p^s^i) depend on the ordinary or ex- 
traordinary refractive indexes Uj of the crystal [18]. The 
TPSA determines all the spectral and temporal proper- 
ties of two-photon light. In the general form, its square 
module gives the spectrum of the biphoton radiation. For 
example, the signal spectrum is 
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duJiF{oJs,uJi 



(3) 



and similarly for the idler spectrum. 

In bulk crystals the momentum conservation is en- 
sured, for each frequency, by choosing a proper direc- 
tion of the waves involved in the process. Nevertheless, 
in type II SPDC the efficiency of the process is reduced 
by birefringence walk-off. The advantage of using poled 
crystals for biphoton generation relies on the absence of 
transverse walk-off and allows one to avoid using narrow- 
band filtering in certain experiments. Moreover, com- 
pared with a bulk crystal, a poled material provides a 
higher SPDC pair production efficiency due to noncritical 
quasi-phase matching (QPM), obtainable at a given tem- 
perature, or due to the possibility to utilize the largest 
value of the effective nonlinear coefficient. 

According to [22| , let us consider a square nonuniform 
first-order QPM grating with slow variation of the poling 
period A{z). Then the spatial variation of the second- 
order susceptibility along the medium can be written as 



(2) iKoz+i^oiz) 

AO ^ 5 



(4) 



where the absolute value of the second-order susceptibil- 
ity is constant. The phase ^o{z) in Eq. dH is induced by 
the modulation of the local wavevector hz{z) = 27r/A(z). 
The wavevector varies linearly within the length of the 
crystal (linearly chirped QPM grating), thus [22| 
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where a is defined as the chirping parameter and Kq is 
chosen so that perfect collinear first-order quasi phase 
matching for selected frequencies ujso and = Up — Uso 
is obtained at the center of the crystal, i.e., 

kp{ujp) - ks{uJso) - kii^p - ^so) - Ko = 0. (6) 



The phase of the grating, obtained by integrating Eq. (j5j) 
(except for an unessential constant phase factor), is 



Mz) = -a{z^^f 



(7) 



and it equals zero at the center of the crystal. 

The TPSA gains a phase shift that depends nonlinearly 
on the frequency because of the presence of this phase 
factor in the second-order nonlinear coefficient. In fact, 
doing analytically the integral in Eq. ([2]), we find that 
the TPSA takes the form 
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xS{ujp - Us - uji) (8) 



where Erf is the error function and Ak{uJs,uJi) = kp — kg — 
ki — Kq. The phase factor shows a nonlinear dependence 
on Ak{uJs^^i)' 



(j){uJs,OJi 



Ak{uJs^uJi)L Ak'^{ujs,uji) 
2 ^ 4^ ' 



(9) 



As a consequence, the phase depends nonlinearly on the 
frequency, as discussed in more detail in section Hill This 
phase factor does not affect the spectrum of the biphoton 
wavepackets as it is evident from (|3|). Nevertheless, it is of 
huge importance when considering the temporal behavior 
of the wavepacket [9] . On the contrary, the error function 
in Eq. ([8j) affects the width of the spectrum of biphotons 
and the effect is the broadening of the spectrum when 
increasing the aperiodicity. 



III. TEMPORAL PROPERTIES 

In Fig. [1] a typical experimental set-up for the mea- 
surement of the time separation between the arrivals of 
the two photons, r, is depicted. Photons of a pair pro- 
duced by type II SPDC are orthogonally polarized, and 
can be easily separated by a polarizing beam-splitter and 
addressed to two different photon counting detectors, Dl 
and D2. 

The rate of coincidences between signal and idler pho- 
tons reaching detectors at times t and t , respectively, 
can be expressed as [13] 

Rc^ JJ dtdt'G^^l,,{t;t')U{t-t\Tw), (10) 



where 11 (t 

dow function of temporal width Tw and Grnias {t ;t j is 
the second-order Glauber's correlation function, where, 
for convenience, not only the electric field operators are 



t jTw) is the rectangular coincidence win- 

(2) 
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FIG. 1: (color online) A typical set-up for observing two- 
photon coincidences in type-II SPDC. A laser beam pumps 
a type-II aperiodically poled crystal (APPC) in collinear 
regime. After a polarizing beam splitter (PBS), signal and 
idler beams are addressed to low-jitter photodetectors, whose 
outputs are sent to a coincidence circuit. 



incorporated, but also the spectral response functions of 
the signal and idler detection arms, fs and f^, respec- 
tively. Gineas ; O stands then for the measured value 
of Glauber's correlation function. Then, 

= |K0|Bi(t")i,(t'^'^'-^"2 

where 



*■'=/ 



A, 



B,{t") = j dujU^,)e-^^^' h^^. 

Typically in the literature (see for example Ref. pjj ) 
there is no spectral dependence of the detection arms, 
and the operators As {t ) and Bs {t ) are usual quantized 
electric fields, i.e.. 



(12) 



with X = However, in the general case the spectral 
response of the detection channel should be taken into 
account since it influences both spectral and temporal 
properties of the measured two-photon wavepackets. For 
instance, considering Gaussian spectral responses of the 
detection channels, we have 



(13) 



where ujxc is the selected central frequency of the spec- 
tral response and 5uj is its width. According to Ref. [23| , 



the measured time two-photon amplitude (TTPA) is ex- 
pressed as 

{0\Bi{t')As{t)\ilj) = J duJsduJiF{ujs.cJi)fs{uJs)fM) 

-^^^^".(14) 
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The measured spectrum (|3]) of the wavepacket is also 
modified by the detection arms spectral filtering as 



S{uJs) 



J duJiF{ujs,oJi)fs{oJs)fi{oJi) 



(15) 



Until the end of the next section we will consider only 
the simplest case of detection arms without any spec- 
tral dependence, i.e., fx{(^x) = 1, in order to focus on 
the temporal properties of the emission. The discussion 
on the modification of the temporal properties induced 
by the spectral characteristics of the detection arms is 
postponed to section |Vl To discuss the difference in the 
arrival times of the signal photon, t , and the idler one, 
t , on the detectors, we introduce the new variables: 



T 

T 



(16) 



Thus, the TTPA can be rewritten as the Fourier trans- 
form of the TPSA: 

F(r) = {0\B,{i)As{tM 

OC J duJsF{uJs,ojp - ujs)e'^^"^ ■ 

The square module of the TTPA is the measured 
Glauber's second-order correlation function. 



(17) 



(2) 

and the width of Grnias 
correlation time. 



(r) will be further referred to as 



IV. TEMPORAL BIPHOTON COMPRESSION 

Since the TTPA is the Fourier transform of TPSA, a 
necessary condition for reaching ultra narrow correlation 
time is producing ultra broadband SPDC light. In the 
case of aperiodically poled crystals with linear chirping, 
the phase factor depending nonlinearly on z is responsi- 
ble for the spectral broadening effect and increasing the 
aperiodicity broadens the spectrum. Unfortunately, the 
two-photon spectral amplitude is not Fourier transform 
limited because of the presence of the phase factor (|9]) in 
Eq. dSj). 

Thus, the spread in the detection time difference of the 
two photons cannot be reduced simply by broadening the 
spectrum. As suggested in Ref. [9] , the ideal way to make 
TTPA perfectly Fourier transform limited is to insert in 



the path of the biphotons a proper optical medium that 
compensates for the nonhnear part of the phase factor in 
Eq. dnj. In general the insertion of an optical medium 
with refractive index rimxi^) and of length in both 
optical paths x = s^i adds to the TPSA the phase factor 
g^(fc^,(a;,)^,+/c^,(a;,)^^)^ whcrc kmx = ""^^^^^^"^^ . Denoting 



the TPSA after the propagation of biphotons through the ^ 



optical medium as F^, we have 



In Ref. [9] it was shown that a correlation time at the 
Fourier transform limit can be achieved when the inserted 
optical media satisfies the condition 

kms{^s)ls + kmi{(^p - (^s)k = -(t){uJs,^p " ^s)- (19) 

In Ref. [l5| it was shown that, in some specific cases, 
this can be achieved by passing one of the two photon 
through a normal-dispersion medium, for example a stan- 
dard optical fiber. This result opens a perspective of a 
real implementation of the method suggested in Ref. @ . 
In particular, we considered the simplest case where in 
Eq. (p!9|) the terms of orders higher than two in the fre- 
quency expansion of the wave vector are negligible. Here 
we explicitly discuss the effect of considering such higher- 
order terms in the frequency expansion of Eq. ([19]). The 
first obvious consideration is that a perfect Fourier trans- 
form Hmited compression does not require Eq. (p!9|) to be 
satisfied; in fact, it is sufficient to satisfy its second-order 
derivative: 



(20) 

Furthermore, it is worth observing that Eq. (|2Q|) should 
be satisfied only for the values of uJs where the TPSA is 
nonzero. 

In order to understand how to satisfy, at least approx- 
imately, the condition in Eq. (|2Q|) we compare the be- 
havior of the second-order derivative of the Harris Phase 
factor, —^^(bi ) (HP), which needs to be elim- 

inated to provide the temporal compression of the bipho- 
ton at the Fourier limit, and the one induced by the pres- 
ence of the optical fibers of appropriate lengths in the sig- 
nal and idler channels, -£j2[kms{^s)ls + ^mi{^p — ^s)h] 
(OFP). We considered collinear type II, almost degen- 
erate, SPDC emission in an aperiodically poled KTP 
crystal (APKTP). The monochromatic pump at 458 nm 
and the idler field are ordinarily polarized, while the sig- 
nal field is extraordinarily polarized. The dependence of 
wavectors on the frequency in KTP is evaluated exploit- 
ing Sellmeier equations [24] and we consider the first- 
order QPM. The phase mismatch at the center of the 
APKTP crystal is compensated for by the inverse grating 
vector Kq = 2441.8 cm~^, corresponding to degenerate 
emission at 916 nm. Spectral and temporal properties of 
the SPDC radiation, as well as the temporal compression 
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FIG. 2: (color online) The dependence of HP and OFP on the 
signal wavelength (both for positive and negative chirping) 
with the optical fiber present only in the signal channel, for a 
crystal with L = 0.8 cm and |a| = 50cm~^. For completeness, 
the spectrum of the biphoton field is shown as well 



of the biphoton wavepackets in the optical fiber, are nu- 
merically calculated for several different parameters such 
as the length of the crystal, the chirping parameter and 
the fiber length. No frequency expansion of the wavevec- 
tor is involved, in contrast to Ref. p^] where the analysis 
was performed analytically and only up to second-order 
terms in frequency. 

Fig. ([2]) shows the dependence of HP and OFP on 
the signal wavelength in the case where the optical fiber 
is present only in the signal channel, for a crystal with 
L = 0.8 cm and |a| = 50cm~^ . For completeness, the 
spectrum of the biphoton field is shown as well. It is in- 
teresting that the effect of the optical fiber in either signal 
or idler SPDC channel is different for the cases of posi- 
tive {a > 0) and negative {a < 0) chirping parameter. In 
fact, since temporal compression is achieved due to the 
group velocity dispersion in a standard optical fiber, only 
in the case of negative chirping the HP can be approxi- 
mately compensated by the OFP. Temporal compression 
in the case of a > can be obtained only in a specially 
engineered fiber with negative group velocity dispersion, 
whose study is beyond the scope of this paper. 

The maximum value of the compression is achieved for 
the length of the fiber // providing the best agreement 
between the HP and OFP, approximately corresponding 
to the situation where HP and OFP have the same value 
at the degenerate wavelength (916 nm), as shown in Fig. 
[21 The dependence of the correlation time on the optical 
fiber length is shown in Fig. [3] both for positive and 
negative chirping parameter. The effect of the optical 
fiber in the first case is an increase in the correlation 

(2') 

time, i.e., broadening of GmeasiT)^ while in the second 
case the insertion of the fiber leads to a compression of 
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FIG. 3: (color online) Dependence of the correlation time 
on the optical fiber length. The two lines refer to positive 
(dashed line) and negative (solid line) chirping parameters. 



FIG. 4: (color online) Cl^iLs (r) evaluated via exact numerical 



Grneas{T), Until its minimum width is reached at / = 
The compression is then followed by a broadening of 

G^^^ (r) 

At / = the computed correlation times do not have 
the same values for negative and positive chirp. This 
effect is due to the generation of signal and idler photons 
at different frequencies at the beginning and at the end 
of the crystal. In fact, biphotons are degenerate only 
at the center of the crystal and the frequency difference 
causes distinct delays between the two photons. During 
the propagation of biphotons in the crystal, the delay 
due to the frequency difference is compensated by the e- 
o delay in the case of a > but in the case of a < 0, the 
delays add up. 

In tl5] it was shown that within the second-order fre- 
quency series expansion of (j){ujs^ujp — cjg), perfect tem- 

poral compression of Grnias {t) can be achieved with nor- 
mal dispersive media. Let us now discuss the effect of 
higher-order terms in the expansion on the width of the 
compressed Grneas{T). Fig. IHreports the results obtained 
via exact numerical calculation for a crystal with L = 0.8 
cm and \a\ =20 cm~^. The left-hand part of the figure 

(2) 

shows the shape of GmeasiT) at the output of the crystal, 
with no dispersive medium inserted, while the right-hand 
part illustrates the effect of a 1.06 m long fiber inserted 
into the SPDC signal channel. The a) and b) parts refer 
to a crystal with a negative chirping parameter, a < 
, while the c) and d) parts to an a > crystal. In the 

case of negative chirping, the compressed Grneas{T) takes 
almost the shape of an ideal sine function, meaning that 
time compression at the Fourier limit is almost achieved 
[ilj. The estimated FWHM of G^^Ls(r) is reduced from 
2.8 ps to 67 fs (which is the value at the Fourier limit), as 
it is clear from Fig. [5l where the spectrum of the corre- 
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^meas V 



calculation for a crystal with L = 0.8 cm and | 
The left-hand part of the figure shows the shape of Gmlas{r) 
at the output of the crystal, with no dispersive medium in- 
serted, while the right-hand part illustrates the effect of a 1.06 
m long fiber inserted into the SPDC signal channel. The a) 
and b) parts refer to a crystal with a < while the c) and d) 
parts, to an a > crystal. 
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FIG. 5: (color online) The spectrum of biphotons, with the 
width of approximately 40 nm, as well as HP and OFP fre- 
quency dependencies (dashed and solid lines, respectively). 



sponding biphoton wavepacket is shown. Summarizing, 
Fig. m shows that a normal dispersive medium, as an 
optical fiber, can compress the biphoton wavepacket to 
the Fourier limit. The effect is opposite in the case of 
a > 0: the initial correlation time of 2.6 ps is increased 
to become 5 ps. 

With this choice of crystal parameters, the condition 
given by Eq. (f2Q|) should be satisfied in a quite narrow 
region of nonzero signal field spectrum. Fig. [5] shows 
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FIG. 6: (color online) The spectrum of biphotons as well as 
HP and OFP frequency dependencies (dashed and solid lines, 
respectively), for a crystal with L = 1.8 cm and a = — 50 
cm~^ 



FIG. 8: (color online) The spectrum of biphotons as well as 
HP and OFP frequency dependencies (dashed and solid lines, 
respectively), for a crystal with L = 2.5 cm and a = — 50 
cm~^ 
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FIG. 7: (color onhne) Effect of inserting a 0.52 m long fiber FIG. 9: (color onhne) Effect of inserting a 0.57 m long fiber 
(right-hand side) on the correlation time for the case shown (right-hand side) on the correlation time for the case shown 
in Fig. dSl) in Fig. [8] 



the spectrum of biphotons, approximately 40 nm wide, 
as well as HP and OFP lines: their superimposition, in 
the case of negative chirping, is easily obtained in such 
a narrow spectral region. However, an increase in either 
the length of the crystal or the chirping parameter leads 
to a broadening of the spectrum. 

Figs. [6] and [8] are analogues of Fig. [5] for two differ- 
ent crystals. In both cases, the same chirping parameter 
OL = —50 cm~^ is considered while the length is L = 1.8 
cm in the first case and L = 2.5 cm in the second one. 
As the crystal length is increased, the spectrum becomes 
broader and broader, its width growing from 300 nm in 



the first case to 500 nm in the second one. As a conse- 
quence, the HP and OFP frequency dependencies must 
be superimposed in a larger region, and, as in general this 
cannot be achieved, the perfect temporal compression of 
the Grneas{T) is sHghtly deteriorated. 

The compressed peaks become more noisy and start to 
differ from the sine function and, as a consequence, the 
width is slightly increased. Figs. [3 and [9] show the effect 
of the insertion of the fiber. The length of the fiber is 
0.52 m in the first case and 0.57 m in the second case. 
However, a significant compression still remains: from 7.5 
ps to 700 fs in the first case and from 11.25 ps to 1 ps in 
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FIG. 10: (color online) Correlation time as a function of 
the length L of the crystal with the chirping parameter 
a — lOOcm"^ (solid line) and a — — lOOcm"^ (dashed line). 

the second case. The correlation time is reduced in both 
cases by about an order of magnitude, confirming that 
ultrabroadband PDC light can be generated in chirped 
crystals and that temporal compression can be achieved 
in a standard optical fiber. 

Fig. [To] shows the correlation time as a function of 
the length L of a crystal with the chirping parameter 
a = lOOcm"^ (solid line) and a = — lOOcm"^ (dashed 
line). Each point refers to the value that the correlation 
time takes after the passage of the signal photon through 
a fiber with the length optimized for the maximum com- 
pression in the case of negative-chirped crystals. The plot 
demonstrates that the difference between the correlation 
times for positive and negative chirping parameters be- 
comes more pronounced for longer crystals. In fact, an 
increase in the crystal length leads to a growth of the 
correlation time in the case of a > 0. Calculation shows 
that for a < 0, an increase in L affects the correlation 
time very little, although a perfect compression would 
result in the correlation time inversely proportional to 
L. In any case, the width of compressed G\neas{T) is al- 
ways below 200 fs, at least for the range of crystal length 
considered. 

It is interesting to point out that the curves in Fig. 
[TOl are not affected by changing the absolute value of 
the chirping parameter. We have computed the corre- 
lation times for both positive and negative a with the 
absolute value varying in the range 20 - 500 cm~'^, for 
three lengths of the crystal (0.8 cm, 1.8 cm and 2.5 cm) 
and no significant difference was observed. This seems 
very promising for the temporal compression of ultra- 
broadband biphotons, as increasing the chirping param- 
eter leads to a broadening of the spectrum, as shown in 
Fig. [TTJ For example, for a crystal of L = 2.5 cm, if 
the absolute value of the chirping parameter is increased 
from 50 cm~^ to 500 cm~^ the corresponding spectrum 



FIG. 11: (color online) Correlation time as a function of the 
chirping parameter, for three different crystal lengths. 

width will change from 53.28 nm to 425.86 nm and, at 
the same time, the correlation time of the compressed 
biphoton wavepacket can be still below 1 ps simply with 
a standard fiber in the signal channel. 

Furthermore, we found that the length of the fiber to 
be inserted for the temporal compression is almost inde- 
pendent of the crystal length. At the same time, it is 
reduced with the increase of a, passing from 2 m to 50 
cm approximately in the range of the chirping parameters 
considered above. 



V. SPECTRAL FILTERING EFFECTS 

In the previous section we considered the spectral and 
temporal properties assuming the absence of spectral se- 
lection in the detection arms. Here, we discuss how the 
presence of such spectral selection associated with the de- 
tection process modifies the measured spectral and tem- 
poral properties of the biphoton wavepacket as well as 
the measurement of its temporal compression. In par- 
ticular, we consider detection arms with Gaussian spec- 
tral filtering as the one in Eq. ([13]) centered at 532 nm 
{ujsc = ^ic = 3222 THz) and width 250 nm {Suj = 716 
THz). Fig. [represents the effective spectral behavior of 
such detection arms for signal/idler wavelengths in the 
region of interest, and it is in agreement with the typi- 
cal spectral response of certain commercial single-photon 
detectors. 

For comparison, we analyze the case of two APKTP 
crystals already considered in the previous section. The 
first one has the chirping parameter a = — 20 cm~^ and 
length L = 0.8 cm. The spectrum not accounting for 
detection spectral filtering is shown in Fig. [5l and the 
corresponding temporal behavior is shown in Fig. [1] The 
second one has parameters a = — 50 cm~^ and L = 2.5 
cm, and its spectral and temporal behavior (with no ac- 
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FIG. 12: (color online) Simulated spectral response of a 
single-photon detector. 
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FIG. 14: (color online) Effect of a 1.06 m long fiber inserted 
into the signal channel (right hand side) on Gm\as (t) for the 
case shown in Fig. [13] 
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FIG. 13: (color online) The modified spectrum of biphotons 
for a crystal with L — 0.8 cm and a — —20 cm~^. 



count for the detection spectral filtering) is shown in Fig. 
[8] and Fig. [9l respectively. Figs. [13] and [15] show the 
spectrum of the biphoton wavepackets modified by the 
detection spectral filtering of Fig. [12] in the case of the 
first and the second crystals, respectively. 

As the spectrum in Fig. [13] is quite narrow and its 
shape is not modified much with respect to the one in Fig. 
ISl it is reasonable to expect that also the corresponding 
Grneas(T), both with and without the fiber for temporal 
compression, is not modified much by the presence of 
the detection spectral filtering. This is confirmed by the 
comparison between Fig. [14] and Fig. |4] 

On the contrary, the wide spectrum of Fig. [8] is 
strongly distorted by the presence of the detection spec- 
tral filtering as shown in Fig. [15] The temporal coun- 
terpart of this distortion can be found in Fig. [16] where 
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FIG. 15: (color online) The modified spectrum of biphotons 
for a crystal with L = 2.5 cm and a = — 50 cm~^. 



Grneas{T) manifests shapes modified with respect to the 
ones in Fig. [9l even if the correlation time widths look 
not so different. 

Summarizing, the spectral filtering in the detection 
arms affects not only the measured spectral properties 
of the biphoton wavepacket, but also its measured tem- 
poral properties. Thus, to achieve temporal compression 
with ultrabroadband PDC light, one has to take care of 
the spectral filtering of the detection systems for its ef- 
fect on both spectral and temporal behavior of biphoton 
wavepackets. 
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FIG. 16: (color online) Effect of inserting a 0.57 m long fiber 
(right-hand side) on Gmlasir) for the case shown in Fig. [15] 

VI. CONCLUSION 

Producing nonclassical light with tailored spectral and 
temporal properties represents a fundamental resource 
for developing quantum technologies, such as quantum 
lithography, quantum metrology, two-photon coherent 
absorption, etc. Aiming at an effective practical real- 
ization of our ideas, first presented in 0, in this paper 
we have systematically discussed the possibility of gener- 
ating two-photon wavepackets with correlation times as 
small as possible by exploiting aperiodically poled crys- 
tals together with standard optical fibers. The analysis 
carried out here is numerical, taking into account the ex- 
act dispersion dependencies in the crystal and the fiber, 
given by the Sellmeier equations. This is a consider- 
able step forward compared to the previous work where 
the effect was mainly described analytically, but only in 
the framework of the second-order approximation in the 
frequency expansion of the phase mismatch. From our 
numerical calculations it follows that higher-order terms 



in the frequency expansion of the wavevector mismatch, 
as well as in the fibre dispersion dependence, consid- 
erably change the shape of the compressed two-photon 
wavepacket. However, the effect of compression is still 
present. The dependence of this effect on the crystal 
length and the chirping parameter has been analyzed. 
Our results clearly show that a significant compression, 
i.e. reduction of correlation times by more than one or- 
der of magnitude, can be easily reached, achieving cor- 
relation times as short as hundreds of femtoseconds. In 
our analysis we also discussed the effect of possible spec- 
tral filtering in detection arms. Despite spectral filtering 
modifies the correlation time of the biphoton wavepacket 
as well as its spectrum, it is noteworthy to observe that, 
at least in the case considered, the temporal compression 
is essentially preserved. Direct experimental observation 
of such a compression by means of coincidence measure- 
ment is possible only if the detectors' jitter is smaller than 
the correlation time of the biphoton wavepacket. Even 
if hundreds of femtoseconds jitter seems hardly achiev- 
able for present days detector technologies, it could be 
a target of next future detectors |25l-[28| . furthermore, 
undirect measurements can be envisaged [29]. Since our 
scheme only requires relatively simple and available tech- 
nologies (i.e. chirped poled crystals and commercial op- 
tical fibers), we have demonstrated the feasibility of a 
suitable practical realization of ultracompressed bipho- 
tons. 
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